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ABSTRACT 

We present the first robust measurement of the high redshift mass-metallicity (MZ) relation at 10 8 < M/Mq < 10 10 , 
obtained by stacking spectra of 83 emission-line galaxies with secure redshifts between 1.3 < z < 2.3. For these 
redshifts, infrared grism spectroscopy with the Hubble Space Telescope Wide Field Camera 3 is sensitive to the 
/?23 metallicity diagnostic: ([On] AA3726, 3729 + [Om] AA4959, 5007)/H/3. Using spectra stacked in four mass 
quartiles, we find a MZ relation that declines significantly with decreasing mass, extending from 12+log(0/H) = 

8.8 at M = 10 9 - 8 M q , to 12+log(0/H)= 8.2 at M — 10 8,2 M o . After correcting for systematic offsets between 
metallicity indicators, we compare our MZ relation to measurements from the stacked spectra of galaxies with 
M > 10 9 - 5 M o and z ~ 2.3. Within the statistical uncertainties, our MZ relation agrees with the z ~ 2.3 result, 
particularly since our somewhat higher metallicities (by around 0.1 dex) are qualitatively consistent with the lower 
mean redshift (z = 1.76) of our sample. For the masses probed by our data, the MZ relation shows a steep 
slope which is suggestive of feedback from energy-driven winds, and a cosmological downsizing evolution where 
high mass galaxies reach the local MZ relation at earlier times. In addition, we show that our sample falls on an 
extrapolation of the star-forming main sequence (the SFR-M* relation) at this redshift. This result indicates that 
grism emission-line selected samples do not have preferentially high star formation rates (SFRs). Finally, we report 
no evidence for evolution of the mass-metallicity-SFR plane; our stack- averaged measurements show excellent 
agreement with the local relation. 
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1. INTRODUCTION 

The correlation between stellar mass and gas-phase metallic- 
ity in galaxies (the mass-metallicity or MZ relation) is sensi- 
tive to the processes that regulate their growth. This relation is 
shaped by gas outflow and accretion rates, the enrichment of 
these gas flows, and the star formation efficiency of galaxies 
(Tremonti et al. 2004; Brooks et al. 2007; Dalcanton 2007; Erb 
2008; Finlator & Dave 2008; Peeples & Shankar 2011; Dave 
et al. 2012; Henry et al. 2013). It is expected, then, that observa- 
tions of the MZ relation toward low masses and high redshifts 
will constrain the assembly of galaxies and the metal enrichment 
of the intergalactic medium. 

To date, a complete picture of metallicity evolution has not 
emerged. While studies focusing on higher mass galaxies at 
0.5 < z < 3 have shown lower metallicities than are seen 
locally (Zahid et al. 2011; Shapley et al. 2005; Erb et al. 2006; 
Liu et al. 2008; Maiolino et al. 2008; Wright et al. 2009; 
Hayashi et al. 2009), evolution below M ~ 10 9 °-10 95 M o 


* Based on observations made with the NASA/ESA Hubble Space Telescope, 
which is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555. 


remains unconstrained. In Henry et al. (2013), we began to probe 
masses just below M ~ 10 8 5 M© at intermediate redshifts; 
nevertheless measuring metallicities at higher redshifts (z > 1) 
has been challenging because of the requirement for infrared 
spectroscopy of large, faint samples of galaxies. Although lower 
masses and metallicities have been reached for some strongly 
lensed galaxies (e.g., Teplitz et al. 2000; Hainline et al. 2009; 
Wuyts et al. 2012; Brammer et al. 2012; Yuan et al. 2013; 
Belli et al. 2013), the uncertainties of these measurements are 
large. 

In this Letter, we present the MZ relation at z > L3, 
derived by stacking 83 galaxies with observations covering 
the R 2 3 metallicity diagnostic ([On] AA3726, 3729 + [Om] 
AA4959, 5007)/H/3 (Pagel et al. 1979). Our sample, which is 
drawn from Hubble Space Telescope Wide Field Camera 3 
(WFC3; MacKenty et al. 2010) grism observations, reaches 
M ~ 10 8 M o at 1 .3 < z < 2.3. By focusing on measure- 
ments from stacked spectra, we avoid bias introduced by the 
requirement that Up be detected in each spectrum. We adopt a 
Chabrier (2003) initial mass function (IMF) and take Q M = 0.3, 
£2 A = 0.7 and H 0 = 70 km s -1 Mpc -1 . When reporting mea- 
surements of doublet lines (i.e., [Om] AA4959, 5007), we use 
“[O ii]” or “[O hi]” to refer to both lines. 
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2. OBSERVATIONS AND SAMPLE SELECTION 

We use data from the WFC3 Infrared Spectroscopic Parallel 
(WISP) Survey, as described by Atek et al. (2010, 2011), 
Dominguez et al. (2013), Colbert et al. (2013), and Bedregal 
et al. (2013). The WISP Survey uses both the G102 (0.8-1 . 1 /zm, 
R ~ 210) and G141 (1.1-1 .7 /zm, R ~ 130) grisms, enabling 
measurements of R 23 from 1.3 < z < 2.3. As the WISP 
survey is ongoing, we restrict our analysis to 29 fields that 
have emission-line catalogs and supporting optical imaging, 
either from WFC3 /UVIS (24 fields) or from the Large Format 
Camera (Simcoe et al. 2000) on the 200-inch Hale Telescope at 
Palomar Observatory ( g and i bands; 5 fields). Five-band galaxy 
spectral energy distributions (SEDs; including Spitzer / IRAC 
3.6 /zm photometry for most) are obtained from photometric 
catalogs described in Dominguez et al. (2013) and Bedregal 
et al. (2013). 

Using the WISP emission line catalogs (Colbert et al. 2013), 
we identify 128 galaxies with secure redshifts measured from 
multiple lines. Objects with contamination from overlapping 
spectra were removed, leaving 74 galaxies. Of these, we exclude 
five galaxies which have inadequate mass constraints: three fall 
in the UVIS chip gap, one is detected in only one broad-band 
filter, and one is blended with another object in ground-based g 
and i images. These cuts leave 69 galaxies from WISP. 

We supplement this sample with publicly released grism 
data covering the Ultra Deep Field (UDF); 12 since the UDF is 
observed only in G141, R 23 is covered only for the 25 galaxies 
between 2.04 < z < 2.36. Of these, we exclude three galaxies 
because of contamination, two because their spectra do not 
contain emission lines, and one with incomplete imaging near 
the edge of the UDF. We construct SEDs for the remaining 19 
galaxies using the nine bands of imaging in the HST extreme 
Deep Field (Illingworth et al. 2013). Photometric catalogs were 
created by using SExtractor (Bertin & Arnouts 1996) in dual 
image mode, with F160W as the detection image. We use 
MAG_AUTO to estimate the total flux. Errors are increased 
to account for correlated noise in drizzled data (Casertano et al. 
2000 ). 

Among the 88 galaxies in the combined sample, we identify 
and remove five active galactic nucleus (AGN) candidates 
(Section 3.3). The mean redshift for the remaining 83 galaxies 
is z = 1.76. 

3. MEASUREMENTS 

3.1. Emission Line Fluxes 

We measure emission line fluxes from the individual spectra 
by simultaneously fitting Gaussian profiles to the [On], [Om], 
and H/3 lines, as well as Ha and [Sn] when they are cov- 
ered. Contribution to Ha from unresolved [Nil] AA6548, 6583 
is taken to be negligible, consistent with our followup spec- 
troscopy (D. Masters et al., in preparation). Each region of the 
spectrum (i.e., [On], [Om] + H/k Ha+ [Sn]) is also fit with a 
linear continuum. Because the dispersion in the grism spectra 
diverges from the adopted solution at the red and blue ends, we 
allow for shifts of up to A z = ±0.01 between the three widely 
separated groups of emission lines. We take the redshift from the 
fit to [O hi] and H/3. Additionally, since the [0 111 ] AA4959, 5007 
doublet is marginally resolved, we fit both lines assuming a 3: 1 
flux ratio. Finally, the emission lines should be spectrally unre- 
solved with widths set by the size of the galaxy in the dispersion 

12 http://monoceros.astro.yale.edu/RELEASE_V3.0/Spectra/UDF/Web/ 
UDF_3dhst_redshift_vL0.html 
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Figure 1 . The MEx diagram shows that a minority of WISP and UDF emission 
line objects are AGN. Filled symbols show the objects where H/3 emission 
is detected at 3cr significance or better (25 galaxies). Open symbols represent 
3cr limits (63 galaxies). The typical stellar mass uncertainty is shown in the 
lower left. Contours show the local relation, derived from SDSS DR7 catalogs 
(http://www.mpa-garching.mpg.de/SDSS). Orange lines show AGN thresholds. 
(A color version of this figure is available in the online journal.) 


direction. Since the dispersion is two times higher in the G102 
spectra, we require lines measured in this grism to have widths 
a factor of two narrower (in A) than those of the redder lines. 

3.2. Stellar Mass 

Stellar masses are derived with the FAST software (Kriek 
et al. 2009), fitting BC03 (Bruzual & Chariot 2003) template 
spectra to our broadband photometry. We take a grid of stellar 
population parameters that include: exponentially declining star 
formation histories with ^-folding times ranging from 40 Myr 
to 10 Gyr; ages ranging from 50 Myr to the age of the universe, 
and A v = 0-3 for a Calzetti et al. (2000) extinction curve. 
Additionally, we use a Chabrier (2003) IMF with metallicities 
of 0.004, 0.008, and 0.02 (Z Q ). We removed the emission 
line contribution to our broadband photometry, estimating the 
corrections by integrating our continuum subtracted spectra (see 
Section 3.4) under the bandpasses that cover the line emission. 

3.3. AGN 

In Figure 1, we use the Mass Excitation (MEx) diagram to 
identify AGN that could contaminate our stacked spectra. This 
diagnostic uses mass as a proxy for the [N 11 ] A6583/Ha ratio, 
which is not measured for our sample. The solid orange lines 
show the AGN threshold from Juneau et al. (2011); objects 
above and to the right are classified as AGN, and objects 
between the solid curves at M > 10 10 M 0 are star- forming/ 
AGN composites. However, the use of the MEx diagnostic is 
compromised by the evolution of the MZ relation. At fixed 
metallicity (or [Oiii]/H/3), galaxies at z ~ 1-2 are 1-2 dex 
more massive than z ~ 0 galaxies (Shapley et al. 2005; Erb 
et al. 2006). Therefore, we expect that the star-forming locus 
in Figure 1 shifts to the right at higher redshifts. The dashed 
curve illustrates the possible evolution, by shifting the Juneau 
et al. curves a conservative 1 dex in stellar mass. Taking this 
modified threshold, we find five AGN candidates. These objects 
are excluded from further analysis. We acknowledge that the 
more aggressive AGN selection would slightly modify our MZ 
relation above M > 10 9,4 M o , increasing the metallicity derived 
in Section 3.5 by 0.03 dex at M = 10 9,8 M 0 . Nevertheless, 
this shift falls within our uncertainties, and does not affect the 
conclusions drawn in Sections 4 and 5. 
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For the galaxies that are undetected in H/3, we cannot rule out 
an AGN for any given object. However, given the distribution of 
detections and limits in Figure 1, it is unlikely that a sizable 
fraction of these galaxies have very high [Oiii]/H/ 3 ratios. 
We also note that none of the UDF galaxies are identified as 
AGN from X-ray counterparts in the Xue et al. (201 1) catalog. 
Therefore, we conclude that stacked spectra will be minimally 
affected by this contamination. Overall, the low AGN content 
in our sample is consistent with the low mass , grism-selected 
objects reported by Trump et al. (2013). 

3.4. Spectral Stacking 

In order to stack spectra, we follow a procedure similar to 
Dominguez et al. (2013). First, we subtract a model continuum 
for each galaxy, made by masking the emission lines and 
smoothing with a 30 pixel boxcar. Each spectrum is then 
normalized by its [Om] flux, so that all galaxies contribute 
equally. Finally, the spectra are de-redshifted and combined by 
taking the median. The relative line fluxes are determined by 
fitting Gaussian profiles. The measurement follows the same 
procedure described in Section 3.1, except the widths of the 
emission lines are allowed to differ from one another because the 
relative contribution of the G102 and G141 data varies with rest- 
wavelength. Error spectra are generated by bootstrap resampling 
the spectra in each mass bin to generate 500 artificial stacks, and 
taking the rms at each wavelength. Errors on the relative line 
fluxes are calculated by propagating the uncertainties on the fit 
parameters. 

Our sample of 83 galaxies is divided into four mass quartiles 
with 7.3 < log M/Mq < 8.6, 8.6 < logM/M 0 < 9.0, 
9.0 < log M/Mq < 9.4, 9.4 ^ logM/M 0 < 10.2. The 
stacked spectra are shown in Figure 2. 

3.5. Metallicity 

A few steps must be taken to infer metallicities from our 
stacked spectra. First, since [Om] and [On] are widely sep- 
arated in wavelength, a dust correction is required. We esti- 
mate extinction by measuring the B aimer decrement for each 
of the five galaxies in our z < 1.5 subsample where both Ha 
and H /3 are detected (at 3 a significance or greater). Consis- 
tent with our findings in Dominguez et al. (2013), we mea- 
sure an average Ha/H/3 = 3.6, with an rms of 1.1. We also 
stack all of the 28 galaxies with 1.3 < z < 1.5; here we find 
Ha/H/3 = 2.48 zb 0.33, consistent with no dust. In the analysis 
that follows, we calculate line ratios assuming both no dust, and 
Ha/H/3 = 3.6 =b 1.1 (with a Calzetti et al. 2000 extinction 
curve). In the latter case, the scatter on the B aimer decrement is 
included in the metallicity errors. 

Second, the flux measured from the H/3 emission line can 
be artificially reduced by stellar absorption, which is estimated 
to have an equivalent width (EW) of 3-4 A in similar grism- 
selected galaxies (Dominguez et al. 2013). Since our stacked 
spectra only represent relative line fluxes, we estimate the H/? 
EW from the mean [O m] EW of the individual stacked galaxies 
and [O hi] /H /3 ratio for that stack. The values for the rest-frame 
H/3 EW range from 30 to 110 A (see Table 1). The relative H/3 
fluxes are increased by (EW res t(H/3) + 3.5)/EW res t to account 
for this absorption. 

Finally, R 23 is double- valued; a single measurement can 
correspond to either a high or low metallicity. Determining 
whether a galaxy lies on the high or low metallicity branch 
can be challenging (Kewley & Ellison 2008; Henry et al. 2013). 



r r : " iL 1 

3500 4000 4500 5000 


Figure 2. Stacked spectra are shown with shaded regions indicating the ±lcr 
uncertainty. Vertical lines mark the expected wavelengths of the line emission, 
including H y + [Oiii]Z4363 and [Nem] A.3869. We do not interpret these 
weaker features because of their low S/N and underlying stellar absorption. 


Nevertheless, under the reasonable assumption that metallicity 
decreases toward lower masses, an average trend is expected: 
as mass decreases, R 23 should increase to log R 23 ~ 1, before 
turning over and decreasing. For the first time at high redshift, 
this trend is observed in Figure 3. We conclude that the 
M = 10 8 8 M q bin marks the “turn-around,” with the two higher 
mass bins falling on the upper branch of R 23 and the lowest 
mass bin implying lower-branch metallicities. The middle panel 
of Figure 3 supports this conclusion; as in the local universe, 
the lowest metallicity galaxies have more highly ionized gas, as 
indicated by higher 032 (log [O iii]/[0 11 ]). 

We have tested that the results from stacking are robust to 
scatter between the upper and lower-branch bins. Taking a 
simulated MZ relation with the same slope as we have measured, 
we adopt an intrinsic metallicity scatter of 0. 1 dex (as at z ~ 0. 1 ; 
Tremonti et al. 2004) and 0.2 dex of scatter in stellar mass due 
to our typical measurement error. We find that fewer than 10% 
of the galaxies in the bins adjacent to the turn-around should 
be attributed to the opposite R 23 branch, and less than 1% of 
galaxies in our high-mass bin will fall on the lower branch. 
This exercise shows that metallicities may be overestimated by 
0.03-0.04 dex in the M = 10 9 2 M© bin, and underestimated 
by the same amount in our lowest mass bin. We do not apply a 
correction for this effect, since its amplitude is smaller than our 
uncertainties. 

Finally, as detailed in Kewley & Ellison (2008), different 
metallicity calibrations of strong-line ratios yield systematically 
different abundances. In Table 1, metallicities are reported on 
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Table 1 

Measurements from Stacked Spectra 


log M*/M© 

N 

(z) 

EW(H/3) 

log [0 111] /H/3 

032 

log R 23 

12+log(0/H) 

12+log(0/H) 

SFR 

12+log(0/H) 




(A) 




(KK04) 

(Maiolino) 

(M© yr- 1 ) 

(FMR) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(ID 

Dust Corrected: Ha/H /3 = 3.6 ± 1.1 

8.2 

21 

1.82 

110 

0.77 ± 0.09 

0.65 ±0.15 

0.85 ± 0.09 

8-16-o‘i5 

7 cq+0.44 
'• Dy -0.16 

7 


8.8 

22 

1.73 

47 

0.81 ±0.08 

0.33 ±0.14 

0.98 ± 0.08 

8.43^ 

8-03! “di 

9 

8.17 

9.2 

21 

1.77 

36 

0.66 ± 0.06 

0.24 ±0.12 

0.85 ± 0.07 

8.65^'-* 


13 

8.33 

9.8 

19 

1.74 

29 

0.53 ± 0.07 

0.22 ±0.13 

0.74 ± 0.08 

0 07+O.O8 

°-° z -o.io 

S.68^ 

21 

8.60 

No Dust Correction 

8.2 

21 

1.82 

110 

0.76 ± 0.09 

0.76 ±0.10 

0.84 ± 0.09 

0 19+O.I5 

12-0.14 

7 C7+0.37 
' —0.15 

3 


0° 

bo 

22 

1.73 

47 

0.82 ±0.08 

0.45 ± 0.08 

0.96 ± 0.08 

8-42^ 0 3 i°6 

8.03t° 0 f 8 

4 

8.22 

9.2 

21 

1.77 

36 

0.67 ± 0.06 

0.35 ± 0.05 

0.83 ± 0.06 

8.70+°,°* 


6 

8.39 

9.8 

19 

1.74 

29 

0.54 ± 0.07 

0.33 ± 0.08 

0.71 ±0.07 

8 - 86 t°o 0 08 

0 77+O.O8 
IZ - 0.09 

9 

8.67 


Notes. (1) Mean stellar mass. (2) Number of galaxies in each stack. (3) Mean redshift in each stack. (4) Average rest-frame H/3 equivalent width of the stacked galaxies 
(uncorrected for stellar absorption). (5) [Oiii]/H/? flux ratio, corrected for stellar absorption. (6) The [Om]/[On] ratio (032 = log [Om]/[On]). (7) The metallicity 
sensitive R 23 diagnostic, corrected for stellar absorption. (8) Oxygen abundance calculated using the KK04 calibration. (9) Oxygen abundance calculated using the 
Maiolino et al. (2008) calibration. (10) SFR estimated from mean H/3 luminosity of the stacked galaxies. (11) The metallicity predicted from the local FMR (Mannucci 
et al. 2011; given for the M > 10 85 M© mass bins where the FMR is defined.) 
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log M/M q 


10.0 


Figure 3. Emission line ratios measured from our stacked spectra (with and 
without dust correction) show that the turn-around metallicity is reached around 
M ~ 10 8 - 8 M o . 

(A color version of this figure is available in the online journal.) 


both the Kobulnicky & Kewley (2004, KK04) and Maiolino 
et al. (2008) scale. 

4. STAR-FORMATION, SAMPLE BIAS, AND THE 
FUNDAMENTAL METALLICITY RELATION 

In the local universe, metallicity, stellar mass, and star forma- 
tion rate (SFR) form a plane called the Fundamental Metallicity 
Relation (FMR; Mannucci et al. 2010), in the sense that galaxies 
with high specific SFRs have lower metallicities. As a result, 
it is important to consider selection effects when interpreting 
the MZ relation. For example, UV-continuum or emission-line 
selected samples may be biased toward higher SFRs and lower 


metallicities than mass- selected samples. Therefore, we use H/3 
luminosity (measured analogously to the H/3 EW, using the 
mean [Om] luminosity and the [Oiii]/H/3 ratio) to estimate 
the SFR for each stack. The H/3 luminosity is corrected for 
dust and stellar absorption, and multiplied by 2.86 to obtain the 
Hoc luminosity. We then calculate the SFR from the relation in 
Kennicutt (1998), converted to a Chabrier (2003) IMF. 13 In 
Table 1 we list the average SFR in each mass bin. Figure 4 shows 
that our dust-corrected SFRs lie very close to an extrapolation of 
the z = L76 star-forming main sequence from Whitaker et al. 
(2012). This comparison shows that our grism- selected sample 
is not likely biased toward high specific SFRs. 

Remarkably, Mannucci et al. (2010) report that the FMR does 
not evolve with redshift (at least out to z ~ 3). Further studies 
confirm this claim for smaller samples of low mass galaxies 
(Henry et al. 2013; Belli et al. 2013). Using our stack- averaged 
measurements, we can also test for evolution. In Table 1 we list 
the oxygen abundance, predicted from the FMR in Mannucci 
et al. (2011). This quantity shows excellent agreement with the 
metallicities calculated on the Maiolino et al. (2008) calibration 
(on which the FMR is defined). We conclude that our sample 
shows no evidence for evolution of the FMR. 

5. THE MAS S-METALLICIT Y RELATION 

Figure 4 shows the MZ relation derived from our stacked 
spectra, for both the dust corrected and uncorrected measure- 
ments. It is clear that the dust correction represents only a small 
adjustment. Excluding our lowest mass bin, our (dust corrected) 
metallicities are 0.33, 0.17, and 0.10 dex lower than the Sloan 
Digital Sky Survey (SDSS) z ~ 0.1 measurement at 10 88 , 10 92 , 
and 10 9 8 M q (Tremonti et al. 2004), suggesting more significant 
evolution in low mass galaxies. A similar trend has been ob- 
served in the intermediate-redshift MZ relation (also shown in 

13 The SFR calibration in Kennicutt (1998) is defined for solar metallicity, and 
will therefore overestimate SFRs at low metallicity (by 25% and 37% for 
continuous star-forming, BC03 models with a Chabrier 2003 IMF at 
Z = 0.4 Z© and Z = 0.2 Z©). However, in order to compare to previous 
measurements of the star-forming main sequence, we report SFRs for solar 
metallicity. 
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Figure 4. Left: the dust-corrected SFRs for our stacked spectra fall along an extrapolation of the z = 1.76 star-forming main sequence, indicated by the dashed line. 
Error bars are derived from the error on the mean Hj8 luminosity, and do not include the uncertainties associated with extinction correction. Right: the MZ relation at 
z ~ 1.76 declines steeply toward low stellar masses. All metallicities are given on the KK04 calibration (when necessary, converted using Kewley & Ellison 2008). 
Intermediate redshift data from Zahid et al. (201 1) and Henry et al. (2013) are binned averages. Orange curves (arbitrarily normalized) show how different prescriptions 
for galactic outflows change the slope of the MZ relation. 

(A color version of this figure is available in the online journal.) 


Figure 4; Zahid et al. 201 1 ; Henry et al. 2013). In fact, in Henry 
et al. (2013) we interpreted this evolution as evidence for down- 
sizing, where higher mass galaxies evolve onto the local MZ 
relation at earlier times. While the uncertainty on the evolution 
measured by the present data remains large, it is consistent with 
downsizing. Stacking larger samples of grism spectra can help to 
reduce the uncertainties, but, at the same time, robust measure- 
ments of the local MZ relation are needed for M < 10 8 5 M©. 
Comparison to local dwarf galaxies remains difficult, as most 
have metallicities derived directly from electron temperatures 
(which differ systematically from strong-line metallicities), or 
suffer from uncertainties regarding the R 23 branch (Berg et al. 
2012; Lyetal. 2013). 

We also compare our MZ relation to the measurement from 
z ~ 2.3 galaxies reported by Erb et al. (2006). These metallic- 
ities were measured from the [Nil] A6583/Ha ratio in stacked 
spectra (using the Pettini & Pagel 2004 calibration), and are 
converted to KK04 metallicities using Kewley & Ellison (2008). 
Figure 4 shows that our MZ relation agrees with the results from 
Erb et al. The somewhat higher metallicities in our data are con- 
sistent with the lower redshift (z ~ 1.76) of our sample. The 
broad agreement between these measurements suggests that the 
systematic offsets between metallicity calibrations (which we 
removed) are not drastically different at low and high redshifts. 
Finally, while some low mass lensed galaxies have metallic- 
ity measurements at these redshifts (e.g., Wuyts et al. 2012), 
large scatter and many weak upper limits preclude meaningful 
comparisons. 

Finally, we compare our data to theoretical models for the 
MZ relation. Following our implementation in Henry et al. 
(2013), we show two models calculated from the prescriptions 
in Dave et al. (2012). Figure 4 shows that a model regulated by 
momentum-driven winds predicts a shallower MZ slope than 
an energy-driven wind model. While we showed in Henry et al. 
(2013) that the intermediate redshift MZ relation is reproduced 
by momentum-driven winds, the higher redshift relation sug- 
gests evolution, as it more closely follows the energy-driven 
wind model. It is unlikely that a different choice of metal- 
licity calibration would flatten our observed MZ relation to 
bring it closer in line with momentum-driven winds; among the 
widely used strong-line metallicity calibrations, only Pilyugin 
(2001) and Pilyugin & Thuan (2005) yield MZ relations 
with significantly shallower slopes compared to a KK04 MZ 


relation. However, these shallow slopes may be inaccurate, ow- 
ing to a saturation of the [O in] A4363 diagnostic at high metal- 
licities (Kewley & Ellison 2008). Nevertheless, the metallic- 
ity errors remain large, so while our data favor energy-driven 
winds, we cannot rule out the shallower slope associated with 
momentum-driven winds. 

We have presented the first measurement of the MZ relation 
from a large sample of low mass, high redshift galaxies. Al- 
though emission-line selected, these galaxies are not extreme 
in their SFRs or metallicities. Rather, our results show excel- 
lent agreement with extrapolations of the star-forming main 
sequence and MZ relation at similar redshifts. We conclude that 
the low-mass, high-redshift MZ relation suggests a downsiz- 
ing evolution and a preference for feedback from energy-driven 
winds. 
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